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Abstract. The article presents a new approach to euclidean plane geometry 
based on projective geometric algebra (PGA). It is designed for anyone with 
an interest in plane geometry, or who wishes to familiarize themselves with 
PGA. After a brief review of PGA, the article focuses on P(M2 ,o,i)i the PGA 
for euclidean plane geometry. It first explores the geometric product involving 
pairs and triples of basic elements (points and lines), establishing a wealth 
of fundamental metric and non-metric properties. It then applies the alge¬ 
bra to a variety of familiar topics in plane euclidean geometry and shows 
that it compares favorably with other approaches in regard to completeness, 
compactness, practicality, and elegance. The seamless integration of euclidean 
and ideal (or “infinite”) elements forms an essential and novel feature of the 
treatment. Numerous figures accompany the text. For readers with the requi¬ 
site mathematical background, a self-contained coordinate-free introduction 
to the algebra is provided in an appendix. 


1. Introduction 

The 19 4/t century witnessed an unprecedented development of geometry and al¬ 
gebra. We need only mention the development of projective and non-euclidean 
geometries, complex and quaternion number systems, and Grassmann algebra to 
indicate the depth and breadth of these developments, many of which came to¬ 
gether in William Clifford’s invention of geometric algebra (jCU78j). This is a 
comprehensive algebraic structure that models both incidence relations and met¬ 
ric relations - for a variety of metric geometries - in a concise and powerful form, 
and which is ideally suited to computational implementation. The teaching and 
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practice of euclidean geometry in the 20 th century, however, remained largely un¬ 
touched by these developments, except for the introduction of vector and linear 
algebra to supplement the standard tools of analytic geometry. 

In recent years, however, geometric algebra has found growing acceptance as 
a tool for euclidean geometry. Those seeking geometric algebra toolkits for doing 
n-dimensional euclidean geometry find two popular solutions in the contempo¬ 
rary literature: the so-called vector geometric algebra (VGA), using n-dimensional 
coordinates (' |DFM07j . Ch. 10); and conformal geometric algebra (CGA), which 
uses (n + 2)-dimensional coordinates (' |DFM07] . Ch. 13). |Gunllb| . [Gunllcj . and 
(Gunlfib] feature a third model, less well known than these two, which fits natu¬ 
rally between them: projective geometric algebra (or PGA for short), which uses 
(n + 1)-dimensional coordinates to model n-dimensional metric spaces of constant 
curvature: euclidean, hyperbolic, and elliptic. This article provides an introduction 
to euclidean PGA, by applying it to the euclidean plane E 2 . 


1.1. Structure of the article 

Sect. [2] gives a brief overview of geometric algebra. Sect. [3] then introduces the dual 
projective geometric algebra P(l &2 o l) as a geometric algebra for doing euclidean 
plane geometry. There follows a discussion of the basis elements in different grades 
and how they can be normalized, along with the distinction between euclidean 
and ideal elements. Sect. [3] examines in detail the geometric product of 2 elements 
of various grades and types, while Sect. [5] does the same for 3-way products. In 
the following sections, the resulting compact and powerful geometric toolkit is 
applied to a sequence of topics in plane geometry: distance formulae (Sect, [d]), 
sums and differences of k -vectors (Sect. [7]), isometries as sandwiches (Sect. 0- or¬ 
thogonal projections (Sect. [9|, a nd a step-by-step solution to a classical geometric 
construction problem (Sect. 10). Sect. [TT] gives the interested reader an overview 
of directions for further study. The article concludes (Sect. 12) by evaluating the 
results obtained and comparing them to alternative approaches to doing euclidean 
plane geometry. Appendix A features a coordinate-free derivation of the results of 
Sect. [3] for readers with the necessary mathematical sophistication. 


2. Geometric algebra fundamentals 

A self-contained introduction to geometric algebra lies outside the scope of this 
article. We sketch here the essential ingredients; interested readers are referred 
to the textbook |DFM07| for a modern computer science approach or Art571 for 
an older, more mathematical approach. The Wikipedia article entitled “geometric 
algebra” is also quite useful. Readers should keep in mind that none of these 
references deal with degenerate metrics, which form a key feature of the approach 
described here. 

Grassmann algebra. Geometric algebra can be built upon the combination of an 
outer and an inner product on a vector space. We assume the reader is familiar 
with real vector spaces, and also with the exterior (or Grassmann) algebra /\ V 
constructed atop a real n-dimensional vector space V. This is a graded algebra in 
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which the elements of grade-fc (/\ fc V) correspond to the weighted vector subspaces 
of V of dimension . Each grade is a vector space in its own right of dimension 
(k). The exterior product 

k m fc+ra 

A: A Vx A v ^ A v 

is a binary operator that is bilinear and anti-symmetric in its arguments. Geometri¬ 
cally, A is the join operator on the subspaces of V: it gives the (fc + m)-dimensional 
subspace spanned by its arguments, or 0 if they are linearly dependent. It is also 
called the outer product. The Grassmann algebra has (n + 1) non-zero grades, 
from 0 (the scalars) to n (the so-called pseudoscalars). /\V has total dimension 
2", as a glance at Pascal’s triangle shows. 

Symmetric bilinear forms. We also assume the reader is familiar with symmetric 
bilinear forms on a vector space, which allow us to define inner products on V. 
Such a form B is characterized by its signature, an integer triple ( p, n, z) where 
p + n + z = n. Sylvester’s Inertia Theorem asserts that there is a basis for V for 
which B is a diagonal matrix with p Vs, n — l’s, and z 0’s on the diagonal. If z ^ 0, 
we say the inner product is degenerate. We will see below that the signature for 
euclidean geometry is degenerate. 

Measurement. In the standard euclidean vector space R 3 , measurement of angles 
between vectors u := (x u , y u , z u ) and v := (, x v ,y v ,z v ) is determined by the stan¬ 
dard euclidean inner product u • v := x u x v + y u y v + z u z v , with signature (3, 0, 0). 
Using this inner product, one can compute the angle between vectors or between 
planes (elements of the dual vector space). If u and v are two unit-length 1-vectors, 
then the inner product u ■ v is well-known to be the cosine of their angle. 

Geometric product. The geometric algebra arises by supplementing the outer prod¬ 
uct with the inner product. One defines the geometric product on 1-vectors of f\ V 
by 

ab := a • b + a A b 

The right-hand side is the sum of a 0-vector (scalar) and 2-vector (plane through 
the origin). This definition can be extended to the whole Grassmann algebra, yield¬ 
ing an associative algebra called the geometric (or Clifford) algebra with signature 
( p,n,z ). In the example above, we obtain R 3i0i o- For details see 1)1 U07 . 

Some terminology. The general element in a geometric algebra is called a multivec¬ 
tor. For a multivector M, the grade-fc part is written (M)*, hence M = ]U fc (M)fc. 
An element of f\ k V is called a fc-vector. A fc-vector that is the product of k 1- 
vectors is called a simple fc-vector, or a blade. For a fc-vector A and an m-vector 
B, the dot product A B := (AB)| fc _ m | is defined as the lowest grade component 
of AB. The wedge A A B = (AB)fc +m is, on the other hand, the highest grade 
component. This is consistent with the definition of the product of two 1-vectors 


uwo elements a and b that satisfy a = Ab for some non-zero A G K represent the same 
subspace, but with different weights. The weight is discussed in more detail below in Sect. 13.4.11 
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Figure 1 . Angles of euclidean lines 

above. X, the reversal of a multi-vector X, is obtained by reversing the ord er o f 
all products involving 1-vectors. X is an algebra involution, needed below in 


5.2 


3. Geometric algebra for the euclidean plane 

The above behavior for R 3 is typical of any geometric algebra with non-degenerate 
metric: the inner product provides the necessary information to calculate the angle 
or distance between the two elements. What is the situation in the euclidean plane 
E 2 ? What kind of inner product do we need to measure the angle between two 
euclidean lines? 

Let 

a 0 x + boy + c 0 = 0, a\X + b\y + c\ = 0 

be two oriented lines which intersect at an angle a. We can assume without loss of 
generality that the coefficients satisfy af + bf = 1. Then it is not difficult to show 
that 

aoai + Mi = cos a 

Unlike the inner product for the case of vectors in R 3 0l o> here the third coordinate 
of the lines makes no difference in the angle calculation: translating a line changes 
only its third coordinate, while leaving the angle between the lines unchanged. 
Refer to Fig. [l] which shows an example involving a general line and a pair of 
horizontal lines. Hence the proper signature for measuring angles in E 2 is (2,0,1). 
This is a so-called degenerate inner product since the last entry in the signature 
is non-zero. 

Notice that to model lines and points in a symmetric way we adopt homoge¬ 
neous coordinates so line equations appear as ax + by + cz = 0. That is, we work 
in projective space KP". Hence, to produce a geometric algebra for the euclidean 
plane we must attach the signature (2,0,1) to a projectivized Grassmann algebra. 
As the above discussion yields a way to measure the angle between lines rather than 
the distance between points, we choose the dual projectivized Grassmann algebra 
P(/\(R 3 )*) for this purpose, where 1-vectors represents lines, 2-vectors represent 
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points, and A is the meet operator. This leads to the geometric algebra P(® 2 ,o 1 ) 
as the correct one for plane euclidean geometry. We call it projective geometric 
algebra (PGA) due to its close connections to projective geometry. (The stan¬ 
dard Grassmann algebra leads to P(R 2 ,o,i)> which models dual euclidean space, a 
different metric space.) 

PGA for euclidean geometry first appeared in the modern literature in [SelOOj 
and |Self)5] and was extended and developed in (Gunllaj , |Gunllb| , [Gunllcj , and 
[Gunl6bj . Readers unfamiliar with duality or projectivization, or just interested 
in a fuller, more rigorous treatment, should consult the latter references. The 
4-dimensional subalgebra consisting of scalars and bivectors, also known as the 
planar quaternions , has a long history as a tool for kinematics in the plane ( )Bla38j . 
[McGDflp . 

3.1. Meet and join 

As mentioned above, the wedge operator A in P(R(; 01 ) ' s the meet operator. It 
is important to have access to the join operator also. Since the typical solution to 
this challenge assumes a non-degenerate metric, we sketch a non-metric approach, 
for details see |Gunllal . The Poincare isomorphism 3 : G G* between the 
Grassmann algebra G and the dual Grassmann algebra G* can be used to define 
the join operator V in P(R 2 o l) : 

A V B := J(J(A) A J(B)) 

J is also sometimes called the dual coordinate map. It is essentially an identity map, 
since it maps a geometric entity in the Grassmann algebra to the same geometric 
entity in the dual Grassmann algebra. For example, in projective 3-space RP 3 , 
a line L can be represented as a bivector in G since it is the join of two points 
(1-vectors in G). It also appears as a bivector in G* since it can also be represented 
as the intersection of two planes (1-vectors in G*). In general, a geometric entity 
represented by a fc-vector in G will be represented by an (n— fc)-vector in G *, where 
n is the dimension of the underlying vector space. J allows one to move back and 
forth between these two dual representations depending on the circumstances. One 
can also implement the join operator using the shuffle operator within P(R)) 0 

(IMIS- Ch. 10). 

3.2. Basis vectors of the algebra 

We provide here a treatment of the algebra based on a choice of basis elements; 
a coordinate-free treatment for more mathematically sophisticated readers can be 
found in Appendix A. 

P(R^ 01 ) has an orthogonal basis of 1-vectors {e 0 ,ei,e 2 } satisfying 
= 0, ef = = 1, e,; • e ? - = 0 for i ^ j 

eo is the ideal line of the plane (sometimes called the “line at infinity”) which we 
write as w, ei is the x = 0 line (vertical!) and e 2 is the y = 0 line (horizontal). All 
lines except u> belong to the euclidean plane and are called euclidean lines. 

We choose the basis 2-vectors 

E 0 := eie 2 , Ei := e 2 e 0 , E 2 := e 0 ei 
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Figure 2. Perspective view of basis 1- and 2-vectors 


for the points of the plane. It is easy to check that these satisfy 
Eg = —1, E 2 = E^ = 0, Ej • Ej = 0 for i ^ j 

Hence the induced inner product on points has signature (0,1,2), more de¬ 
generate than that for lines. As a result, the distance function between points 
cannot be obtained from the inner product - but can be obtained via the geo¬ 
metric product; see Sect. |4.3| below for details. Points that lie on w are said to be 
ideal. Then E 0 is the origin of the coordinate system, Ei is the ideal point in the 
^-direction and E 2 is the ideal point in the y-direction. In general, ideal elements 
can be characterized as elements satisfying x 2 = 0. See Fig. [ 2 ] for a perspective 
view of the fundamental triangle determined by these elements. 

The basis vectors chosen above assume that the first coordinate is the homo¬ 
geneous coordinate. This assumption is helpful when stating results that should be 
valid for general dimensions. On the other hand, existing usage often follows the 
opposite convention; for example, the line with equation ax + by + cz = 0 appears 
in the algebra as ce 0 + ae 1 + be 2 . When writing elements of the algebra as tuples, 
we take into account this existing usage. We write the 1-vector m = ceo + aei -\-be 2 
as [o, 6, c] (square brackets), and the 2-vector P = ccEi + yE 2 + ^E 0 as ( x,y,z ) 
(standard parentheses). 

The pseudoscalar I := eoeie 2 generates the grade-3 vectors. It satisfies I 2 = 
0. This is, the inner product, or metric, is degenerate. A 3-vector p has the form 
al for 06 R. While in a non-degenerate metric the magnitude a of a pseudoscalar 
p can be obtained, up to sign, as pi, this is not possible with a degenerate metric 
(since pi = al 2 = 0 for all p), and we define the signed magnitude S(p) := a. 
We occasionally use the fact for a 1-vector a and a 2-vector P, 5(aAP) = aVP. 
(This follows from the fact if x A y is a pseudoscalar, x V y is a scalar with the 
same magnitude.) Note that Ej was chosen so that ejEj = I. 
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3.3. The geometric product 

The full multiplication table for the basis elements of P(R 2 o 1 ) can be found in 
Table [l] The presence of 0’s indicates that the metric is degenerate, ft is useful to 
have special symbols for the different grade components of the product of two 
blades, which we now provide. Let A be a /c-vector and B, an m- vector. All 
combinations of (k,m) in P(R 2 o 1 ) except (2,2) can then be written as 

AB = A B + A A B 

For (k,m) = (2,2), (AB) 2 =: A x B (= AB-BA), sometimes called the commu¬ 
tator or cross product. We’ll see below that A x B is the ideal point perpendicular 
to the direction of the joining line of A and B. 

3.4. Normalized points and lines 

A k- vector whose square is ±1 is said to be normalized. Since normalization sim¬ 
plifies the subsequent discussion, we introduce it here, although logically speaking 
the justification for all the steps in the normalization process will only later be 
established. The square of any A;-vector in the algebra is a scalar, since all fc-vectors 
in this algebra are simple. Squaring this product and rearranging terms, one ob¬ 
tains a product of the squares of these 1-vectors, each of which reduces to a scalar. 
For a euclidean line m = ceo + aei + 6e 2 , define the norm 

|jm|| := V m 2 = \/m ■ rn (= \J a 2 + b 2 ) 

Then m„ := ||m|| -1 m satisfies m 2 = 1. For a euclidean point P = zEo + xEi + 
yE 2 , P 2 = — z 2 . Define ||P|| := z. Note that, in contrast to a standard norm of 
a vector space, ||P|| can take on positive and negative values, a feature that is 
occasionally useful. Then P„ := z _1 P satisfies ||P„|| = 1. Such a point is also 
called dehomogenized since its Eo coordinate is 1. Note that we have shown that 
normalized euclidean lines have square 1 while normalized euclidean points have 
square -1. In the following discussions we often assume that euclidean lines and 
points are normalized. 
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Table 1. Geometric product in P(R 2 o l) 
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Grade 

Coord. & tuple form 

Norms 

Domain 

Description 

1 

m = aei + be 2 + ceo 
[a, 6, c] 

m | := \Ja 2 + b 2 

IHI ^0 

Euc. line 

llmlloo := c 

IHI = 0 

Ideal line 

2 

P = xEi + j/E 2 + 2 E 0 
(x,y,z) 

IIPII 

IIPIloo 

Euc. point 

IIPIloo := \A 2 + y 2 

l|P||=0 

Ideal point 

3 

al 

11 11 OO •- & 

-all- 

Pseudoscalar 


Table 2. Coordinate-based overview of the euclidean and ideal 
elements with their corresponding norms. 


3.4.1. Weight and norm. If one has chosen a standard representative X for a pro¬ 
jective fc-vector, and Y = AX, we say that Y has weight X. We usually choose 
the standard element to have norm ±1. Such elements of weight ±1 are exactly 
the normalized elements discussed above. The weight can be any non-zero real 
number; while the norm is sometimes restricted to take non-negative values (see 
Table [2] below). The freedom to choose the weight is a consequence of working in 
projective space, since non-zero multiples of an element are all projectively equiv¬ 
alent. Sometimes the weight is irrelevant, sometimes crucial. When multiplying 
elements together, one gets the same projective result regardless of the weights; 
while adding elements, different weights give different projective results. 

3.4.2. Ideal elements and free vectors. Ideal points correspond to euclidean “free 
vectors” (a fact already recognized in |Cli73| ). Let P = aEi + &E 2 be an ideal 
point. Then, as noted above, ||P|| = 0. This leads us to introduce a second norm 
for ideal points, one that is compatible with their function as free vectors. Define 
the ideal norm 

IIPIloo := IIP V Q|| 

where Q is any normalized euclidean point. Then a direct calculation yields 
IIPIloo = \Ja 2 + b 2 , as desired. Thus, the points of the ideal line can be treated as 
free vectors with the positive definite inner product of R 2 (signature (2,0,0)). 

We write the corresponding inner product between two ideal points U and 
V as (U, V)oo. Every euclidean line m has an ideal point moo, normalized so 
that 11 moo 11 00 = 1- The ideal norm allows us to represent ideal points in the 
accompanying figures as familiar free vectors (arrows labeled with capital letters), 
see Fig. [8] (right). 

We also define an ideal norm for ideal lines (i. e., lines m satisfying m 2 = 0). 
For m = aei + t>e 2 + ceo, ||m||oo = c. (As with ||P|| above, this can also take 
on positive and negative values.) Then m = ojj. c > 0 corresponds to an ideal 
line in clockwise orientation; c < 0, to counter-clockwise orientation. Finally for 
completeness we can also consider the pseudoscalar signed magnitude S'(al) as 
an ideal norm: HalHoo := S'(al) = a. We have thus defined an ideal norm for all 
ideal elements in the algebra. This ideal norm, restricted to the ideal plane, has 
signature (2,0,0); considered projectively, this is an elliptic line P(K 2 0 0 ), while 
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Figure 3. Selected geometric products of blades. 


considered as a vector space, it is R^o.o, the geometric algebra of R 2 . See Table [ 2 ] 
for an overview of the euclidean and ideal elements and norms with their domains 
of validity. 

In the following, we will more than once confirm that the standard and ideal 
norms form an organic whole. For a fuller discussion of the ideal norm see §4.4.4 
of [Gunllaj . 

Whether to apply the standard or ideal inner product presents no difficulties 
for practical implementation, as a point can be easily identified as ideal by the 
linear condition P A ui = 0. There is also little danger that an ideal point will be 
mistaken for a euclidean point - all the computational paths that produce ideal 


points presented in this article (see for example Sect. 4.1 Sect. 4.2.2 Sect. 4.3.1 


and Sect. 7.2) produce exact ideal points. This situation is analogous to traditional 


vector algebra: one has no trouble distinguishing vectors and points. 


4. The geometric product in detail: 2-way products 

In the following discussion, P and Q are normalized points (either euclidean or 
ideal, as indicated), and m and n are normalized lines. We analyze the geometric 
meaning of products of pairs and triples of fc-vectors of various grades, paying 
particular attention to the distinction of euclidean and ideal elements. A selection 
of these products is illustrated in Fig. [3] 

4.1. Product with pseudoscalar 

First notice that the pseudoscalar I commutes with everything in the algebra. For 
a euclidean line a, the polar point a 2 - := al = la is the ideal point perpendicular to 
the line a. We can use the polar point to define a consistent orientation on euclidean 
lines; we draw the arrow on an oriented line m so that rotating it by 90° in the 
CCW direction produces nr 1 . See Fig. [ 2 J which shows the resulting orientations on 
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Figure 4. Geometric product ab of two intersecting lines (left) 
and two parallel lines (right). 


the basis 1-vectors. When a is normalized, so is a -1 , another confirmation that the 
two norms (euclidean and ideal) have been harmoniously chosen. For a normalized 
euclidean point P, P x := PI = IP = —e 0 , the ideal line with GW orientation. 
The polar of an ideal point or line is 0. 

We noted above in Sect. 13.21 that the condition I 2 = 0 means the metric is 
degenerate, or, what is the same, multiplication by I (the so-called metric polarity) 
is not an algebra isomorphism. Although some researchers see this as a flaw in 
the algebra (for example, |Li08| . p. 11), our experience leads to view it as an 
advantage, since it accurately mirrors the metric relationships in the euclidean 
plane. For example, when m and n are parallel, = n 2 -, that is, parallel lines 
have the same polar point. In a non-degenerate metric, however, different lines 
have different polar points. In contrast, the degenerate metric properly mirrors 
this euclidean phenomenon. For a fuller discussion of this theme, see Sec. 5.3 of 
|Gunl6bj . 

4.2. Product of two lines 

In general we have mn = (mn)o + (mn )2 = m • n + m A n. We say two lines are 
perpendicular if m • n = 0 - even when one of the lines is ideal. The meaning of 
the two terms on the right-hand side depends on the configuration of m and n as 
follows. 

4.2.1. Intersecting euclidean lines. We say that two intersecting euclidean lines 
meet at an angle a when a rotation of a around their common point brings the 
first oriented line onto the second, respecting the orientation. Then m • n = cos a 
and m A n = (sina)P where P is their normalized intersection point. Consult 
Fig- m left. Readers who are surprised that the angle a can be deduced from the 
wedge product - which doesn’t depend on the metric - are reminded that this is 
possible only because we have used the inner product to normalize the arguments 
in advance. Without normalizing m and n, the formulae are 

m • n = ||m|j|jn|j cosa and m A n = ||m|||jn||(sina)P 

Similar extensions involving non-normalized arguments could be made for the 
subsequent formulae given below, but in the interests of space we omit them. 
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Figure 5. Left: product PQ of two euclidean points; Right: prod¬ 
uct aP of euclidean line and point. 


Exercise: (mn) n = cos not + (sinna)P. Show that the vector subspace generated 
by 1 and P is isomorphic to the complex plane C. 

4.2.2. Parallel euclidean lines, m • n = ±1. We say the lines are parallel when this 
inner product equals 1, otherwise we say they are anti-parallel. In the latter case, 
replace n by —n to obtain parallel lines. Then m ■ n = 1 and m A n = dmnnioo, 
where d mn is the oriented euclidean distance between the two lines and moo- See 
Fig. 0 right. The simplicity of this formula validates the choice of the norm HHoo 
on ideal points. Note that the geometric product in PGA automatically finds the 
correct form of measuring the “distance” between the two lines: the weight of the 
intersection point m A n reflects angle measurement (sin a) for intersecting lines 
and euclidean distance measurement (d mn ) for two parallel lines. 

Exercise: (mn) n = 1 + ndmuUioo. 

4.2.3. Product of a euclidean line with the ideal line. Let n = u be the ideal line. 
Then m n = 0 and mAn = moo is the ideal point of m. Note that since m-n = 0, 
the ideal line is perpendicular to every euclidean line; since it shares an ideal point 
with each such line, it is parallel to every euclidean line! 

4.3. Product of two points 

Here the general formula is PQ = (PQ) 0 + (PQ )2 = P-Q + Px Q. The resulting 
behavior is characterized by the fact that the inner product for points is more 
degenerate than that for lines. 

4.3.1. Two euclidean points. P • Q = —1 and P x Q is an ideal point perpendicular 
to P V Q. To be exact P x Q = — (P V Q)I (notice the negative sign). We also 
write this as (P-Q) x since the ideal point P — Q, rotated in the CCW direction 
by 90°, yields P x Q. See Fig. [5] left. 

Exercise: The distance dpq between two euclidean points satisfies 

dp Q = ||PxQ|! (= ||P V Q||) 
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4.3.2. Euclidean point and ideal point. If Q is ideal, then P • Q = 0 and P x Q 

is the ideal point obtained by rotating Q 90° in the CW direction. This result is 
consistent with the characterization of the product of two euclidean points: it is 
an ideal point perpendicular to PV Q. Q x P rotates in the CCW direction. Thus, 
multiplication of an ideal point by any finite point rotates the ideal point by 90°; 
the specific location of the euclidean point plays no role. 


4.3.3. Two ideal points. The product of two ideal points is zero. Hence the only 
interesting binary operation on ideal points is addition. In light of Sect. 3.4. 1[ this 
helps to explain why ideal points are often treated as vectors rather than projective 
points. 


4.4. Product of a line and a point 

The general formula is mP = (mP)! + (mP) 3 = m • P + m A P. The wedge 
vanishes if and only if P and m are incident. As before, we assume that both m 
and P are normalized. 


4.4.1. Euclidean line and euclidean point, m • P is the line passing through P 
perpendicular to m (consult Fig. [5} right). Why? This can be visualized as starting 
with all the lines through P and removing all traces of the line parallel to m, leaving 
the line perpendicular to m. It has the same norm as m, and its orientation is 
obtained from that of m by CCW rotation of 90°. This is reversed in the product 
P • m. This sub-product is important enough to deserve its own symbol. We define 

nip := m P = —P ■ m 

The wedge product satisfies m A P = d m pl, where d m p is the directed distance 
between m and P. 


4.4.2. Euclidean line and ideal point. Let a be the angle between the direction of 
m and P: cosa = (niocP)^. Then m-P = (cosa)w and mAP = (sina)I. Notice 
that mP is the sum of an ideal line and a pseudoscalar: no euclidean point or line 
appears in the product. The first term, involving the ideal line, is non-zero when 
the ideal line is the only line through P perpendicular to m. When a = every 
line through P is perpendicular to m, and m • P = 0 while m A P = I. 


5. The geometric product in detail: 3-way products 

Products of more than 2 fc-vectors can be understood by multiplying the factors 
out, one pair at a time. The product of 3 different euclidean points (or lines) 
is important enough in its own right to merit a separate discussion. The results 
provide a promising basis for a future investigation of euclidean triangles. Later we 
will see that euclidean reflections (Sect. |8J| ) and orthographic projection (Sect. [9]) 
can also be understood as 3-way products in which one of the factors is repeated. 
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Figure 6. Products of 3 euclidean points 


5.1. Product of 3 euclidean points 

Let the 3 points be A, B, and C. See Fig. [ 6 ] Then using the results obtained above 
for products of two points: 


ABC = (AB)C 

= (-l + (A-B) ± )C 
= —C — (A — B) 

=A-B+C 

The first and second steps follow from the results from Sect. |4.3[ The final equation 
indicates the projective equivalence of the two expressions, since multiplying by 
— 1 does not effect the projective point. The result is somewhat surprising, since 
the scalar part vanishes. Hence, if one begins with the triangle ABC and generates 
a lattice of congruent triangles by translating the triangle along its sides, then the 
vertices of this lattice can be labeled by products of odd numbers of the vertices 
A, B, and C (Fig. [g]). 

Exercise: The product of an odd number of euclidean points is a euclidean point 
that is the alternating sum of the arguments. 


5.2. Product of 3 euclidean lines 

Let the 3 (normalized) lines be a, b, and c oriented cyclically. See Fig. [7] These 
three lines determine a triangle. Then a A b = sin( 7 r —y)C, etc., produces the 
interior angle 7 and the (normalized) vertex C of the triangle. Using the results 
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obtained above for products of two lines: 


abc = (ab)c 

= (— cos 7 )c + (sin 7 )(Cc) 

= (— cos 7 )c + (sin 7 )(C ■ c + C A c) 
= — ((cos 7 )c + (sin 7 ) 0 ^,) + sin 7 d C cI 


The first step follows from the results from Sect. |4.2[ the second and third from 
Sect. 


4.4 


Let C be the intersection of c and C c. In the last equation the expression 
in parentheses is the grade-1 part of the product: b := (abc)i. It is, by inspection, 
minus the result of rotating c around C by 7 . Parenthesizing in a different order 
yields: 

abc = a(bc) = —((cosa)a— (sina)a^) +sinadA a I 


In this form, b is minus the result of rotating a around A by —a. Hence b must 
be the joining line of A and C. See Fig. [tJ 

Since the grade-3 parts are equal, one obtains: 

(sinyjdcc = (sina)d A a 

This illustrates an important technique for generating formulas in geometric alge¬ 
bra. By applying the associative principle one can insert parentheses at different 
positions: 

(ab)c = abc = a(bc) 

The left-hand side and right-hand side represent different paths in the algebra to 
the same result, and these often produce non-trivial identities as this one. 



Figure 7. Product of 3 euclidean lines. 
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Exercises. 1) (abc)i = |(abc + cba). 2) i(cab + cba) = cos( 7 )c. 3) Define 

s := abc + acb -f bac + bca + cab + cba 

Show that s is a 1-vector, called the symmetric line of the triple {a, b, c}. Fig. [7] 
illustrates these relations, and illustrates how the geometric product in PGA pro¬ 
duces compact and elegant expressions for familiar triangle constructions. 


6. Distance and angle formulae 

We collect here the various distance formulae encountered in the process of dis¬ 
cussing the 2-way vector products above. P and Q are normalized euclidean points, 
U and V are normalized ideal points, and m and n are normalized euclidean lines. 
Space limitations prevent further differentiation with respect to signed versus un¬ 
signed distances. Consult Fig. [3] 

1. Intersecting lines. Z(m, n) = cos -1 (m • n) = sin^ 1 (||m A n||) 

2. Parallel lines. d(m, n) = ||m A n||oo 

3. Euclidean points. d( P, Q) = IIP V QII = IIP x Qll^ 

4. Ideal points. Z(U, V) = cos" 1 ((U, V)^) 

5. Euclidean line, euclidean point. d(m, P) = — cZ(P, m) = S(m A P) = m V P 

6. Euclidean line, ideal point. Z(m, U) = cos -1 (||m ■ U||oo) 

Notice that a single expression in the geometric algebra produces several correct 
variants which take into account whether one or the other or both of the arguments 
are ideal. For example, ||m A n|| produces the intersection point of the two lines 
weighted by either the inverse of the sine of the angle (when the lines intersect), or 
the euclidean distance between them (when they are parallel). Similar phenomena 
reveal themselves also in the next section. 

7. Sums and differences of points and of lines 

Based on the discussion of the geometric product above, it is instructive to ex¬ 
amine sums and differences of points, resp. lines. This deceptively simple theme 
reveals important distinctions between euclidean and ideal points and lines that 
play a central role throughout this algebra. It also highlights how traditional vector 
algebra can be directly accessed within P(I &2 0 ,i) ( as the weighted ideal points). As 
before, all points and lines are assumed to be normalized unless otherwise stated. 
Consult Fig. [8] 

7.1. Sums and differences of lines 

When m and n are both euclidean, and intersect in a euclidean point, then m + n 
is their mid-line, the line through their common point m A n that bisects the angle 
between m and n. m — n also passes through their common point, but bisects the 
supplementary angle between the two lines. (To establish the claim, consider the 
inner product of m ± n with each line separately.) If the two lines are parallel, 
then m + n is their mid-line: the line parallel to both, halfway in between them, 
m — n is the ideal line, weighted by the signed distance between the lines. If m 
is euclidean and n = Xuj is a weighted ideal line, then m + n is a (normalized) 
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euclidean line representing the translation of the line m by a signed distance A 
in the direction perpendicular to its own direction (to be exact, in the direction 
opposite its polar point m x ). 

7.2. Sums and differences of points 

When P and Q are both euclidean, P + Q is their mid-point. ( PxC * is the normal¬ 
ized mid-point.) P — Q is an ideal point representing their vector difference. If P 
is normalized euclidean and V is ideal (not necessarily normalized), then P ± V is 
a (normalized) euclidean point representing the translation of the point P by the 
free vector ±V. If both U and V are ideal (again, not necessarily normalized), 
then U ± V is the ideal point representing their vector sum (difference). Here we 
once again meet the R 2 vector space structure on the ideal line induced by the 
ideal norm. 

8. Isometries 

Equipped with our detailed knowledge of 2-way products we now turn to discuss 
how to implement euclidean isometries in the algebra. Recall that the group of 
isometries of E 2 is generated by reflections in euclidean lines. The product of an 
even number of reflections yields a direct (orientation-preserving) isometry (either 
a rotation or a translation), while an odd number produces an indirect (orientation- 
reversing) isometry. Also recall, that in the euclidean plane, every isometry can be 
written using 1, 2, or 3 reflections. We now show how to implement reflections using 
the geometric product, then extend this result to products of 2 and 3 reflections. 

8.1. Reflections 

Suppose a and b are two normalized euclidean lines, and let R a (h) represent the 
reflection of b in a. Purely geometric considerations imply that I? a (b) is a line x 
satisfying a • x = a • b and a A x = b A a. 

Exercise: Show that x := aba fulfils both conditions, satisfies x/b when a/b 
and hence is the desired reflection. 



Figure 8. Left: Sums and differences of normalized euclidean 
lines. Right: Sums and differences involving ideal points and nor¬ 
malized euclidean points. P + Q is the (non-normalized) midpoint 
of segment PQ; P x Q is the ideal point P Q rotated 90° CCW. 
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Figure 9. The reflection in the line a is implemented by the 
sandwich aXa; the product of the reflection in line a followed 
by reflection in (non-parallel) line b is a rotation around their 
common point a A b through 2 cos -1 (a • b). 


Notice that a reflection can then be seen as a special form of a 3-way product 
in which the first and third term is the same line. We write the reflection operator 
b —► aba as a(b). We sometimes refer to this as a sandwich operator since the a 
“sandwiches” the operand b on both sides. 

Exercise: Show that a(P) is also a reflection applied to a euclidean point P. [Hint: 
Write P = mn for orthogonal m and n.] 

8.2. Product of two reflections 

Before we discuss the product of several reflections, we introduce some terminology. 
The product of any number of euclidean lines is called a versor; the product of 
an even number is called a rotor. Versors and rotors are important since sandwich 
operators based on them yield euclidean isometries. 

The concatenation of two reflections in lines a and b can be written 

b(a(x)) = b(axa)b = (ba)x(ab) 

where the expression on the right is obtained by applying associativity to the 
middle expression. Define r := ba, and an operator r(x) := rxr which repre¬ 
sents the composition of these two reflections expressed using the rotor r. Such a 
composition can take two forms, depending on the position of the lines. 

When the lines intersect in a euclidean point, then r is a rotation around 
that point by twice the angle between the lines. See Fig. [9] When the lines are a 
parallel, r is a translation by twice the distance between the lines in the direction 
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Figure 10. Glide reflection generated by r = m + AI applied to 
line x. 


perpendicular to the direction of the lines. The details can be confirmed by ap¬ 
plying the results above involving products of two lines in Sect. T2 to write out r 
for these two cases and then by multiplying out the resulting sandwich operators. 
The rotor for a rotation is called a rotator ; for a translation, a translator. 
Exercises: 1) Show that for a translator t, tx = xt represents half the translation 
of the sandwich t(x). That is, translators also make good “open-faced” sandwiches. 
2) Discuss the rotator cos a + (sina)P when a = |. 


8.3. Product of 3 reflections 

First, recall that a glide reflection is an isometry formed by a reflection in a eu¬ 
clidean line (the axis of the glide reflection) and a translation parallel to this line 
(the order of execution doesn’t matter, since the two operations commute). We 
begin by showing that the sandwich operator generated by the sum of a 1-vector 
and a 3-vector (line and pseudoscalar) corresponds to a glide reflection along the 
fine. Let r = (r)i + (r )3 = m + AI where m is normalized. Then for an arbitrary 
fine x: 

r(x) = rxr 

= (m + AI)x(m — AI) 

= mxm + mxAI — Alxm — A 2 I 2 
= mxm + Amx 1 — Ax x m 
= m(x) + A(mx 2 — x^m) 

= m(x) + 2A(m ■ x x ) 

= m(x) + 2A(cosa)w 

The steps in the calculation follow from the discussion of the 2-way products 
above. The result consists of two terms. The first term is the reflection of x in the 
fine m; by Sect. [7] above, the second term represents the translation of the reflected 
fine perpendicular to its own direction by the distance 2Acos(a). The translation 
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component reveals itself more clearly by considering r(X) for an arbitrary point 
X. A calculation similar to the above yields: 

r(X) = ... = m(X) + 2A(m A X x ) 

= m(X) + 2A(m A cu) 

= m(X) + 2A(m 00 ) 


In this form it is clear that the translation component is 2Am 00 : a translation in 


the direction of the line m through a distance 2A. Consult Fig. 10 


Applying this to the situation of 3 reflections: By Sect. |5.2| above, the product 
of three lines has the form r = abc = b + sin (a)d a Al, hence the above results can 
be applied. Recall that b is the joining line of A and C, the feet of the altitudes 
from A and C, resp. Refer to Fig. [Tj 


8.4. Exponential form for direct isometries 

It’s not necessary to write a rotator as the product of two lines. If one knows the 
desired angle of rotation, one can generate the rotor directly from the fixed point 
P of the rotation. We know that it is normalized so that P 2 = — 1. Then, using 
a well-known technique of geometric algebra, one looks at the exponential power 
series e <p and shows, in analogy to the case of complex number i 2 = — 1, that 
e tp = cost + (sint)P. The right-hand side we already met above as the product of 
two euclidean lines meeting in the point P at the angle t. Setting t = a one obtains 
the rotor r from the previous paragraphs. What’s more, letting t take values from 0 
to a one obtains a smooth interpolation between the identity map and the desired 
rotation. Note that this sandwich operator rotates through the angle 2a; to obtain 

a p 

a rotation of a around P, set r = e~^~ . 

Exercise: Carry out the same analysis for an ideal point V to obtain an exponential 
form for a translator that moves a distance d in the direction perpendicular (CCW) 
to V. [Answer: = 1 + f V.] 


9. Orthogonal projections and rejections 

When one has two geometric entities it is often useful to be able to express one in 
terms of the other. Orthogonal projection is one method to obtain such a decom¬ 
position. For example, in the familiar euclidean VGA any vector b can be 

decomposed with respect to a second vector aasb = aa+/3a x where a, /? £ K. and 
a x ■ a = 0. These two terms are sometimes called the projection , resp., rejection of 
b with respect to a. The algebra P(R 2 o l) offers a variety of such decompositions 
which we now discuss, both for their utility as well as to gain practice in using the 
geometric product introduced above. We can project a line onto a line or a point; 
and a point onto a line or a point. As before all points and lines are assumed to 
be normalized. Consult Fig. El 

Each projection follows the same pattern: take a product of the form XYY 
and apply associativity to obtain X(YY) = (XY)Y. Assuming normalized argu¬ 
ments, YY = ±1, yielding X = ±(XY)Y. The right-hand side typically consists 
of two terms representing an orthogonal decomposition of the left-hand side. Note 
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Figure 11. Orthogonal projections (1. to r.): line m onto line n, 
line m onto point P, and point P onto line m. 


that, like the reflection in a line (in which the first and last factors are identical), 
such projections can be considered as a special form of a 3-way product, in which 
either the first two or the last two factors are identical. 


9.1. Orthogonal projection of a line onto a line 

Assume both lines are euclidean. Multiply the equation mn = m ■ n + m A n with 
n on the right and use n 2 = 1 to obtain 

m = (m • n)n + (m A n)n 

= (cosa)n+ (sinQ-)Pn 
= (cosa)n — (sinci')np 

Note that Pn = -np since P A n = 0. Thus one obtains a decomposition of m as 
left. 

Exercise: If the lines are parallel one obtains m = n + d mn ui. 


the linear combination of n and the perpendicular line np through P. See Fig. 
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9.2. Orthogonal projection of a line onto a point 

Assume both point and line are euclidean. Multiply the equation mP = m • P + 
m A P with P on the right and use P 2 = — 1 to obtain 

m = — (m • P)P — (m A P)P 

= —mpP - (d mP I)P 
— m p diyip co 


In the third equation, nip is the line through P parallel to m, with the same 
orientation. Thus one obtains a decomposition of m as the sum of a line through 
P parallel to m and a multiple of the ideal line (adding which, as noted above in 
Sect. [7} translates euclidean lines parallel to themselves). See Fig. 11 middle. 


9.3. Orthogonal projection of a point onto a line 

Assume both point and line are euclidean. Multiply the equation 


mP = m P + m A P 
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on the left with m on the right and use m 2 = 1 to obtain 

P = m(m • P) + (m A P) 

= m(mp) + m(ci mP I) 

= Pm + dmPm 1 

= Pm + (P _ Pm) 

In the third equation, P m is the point of m closest to m. The second term of the 
third equation is a vector perpendicular to m whose length is dp m : exactly the 
vector P — P m . Thus one obtains a decomposition of P as the point on m closest 
to P plus a vector perpendicular to m. See Fig. El right. 

Exercise: Show that the orthogonal projection of a euclidean point P onto another 
euclidean point Q yields P = Q + (P — Q). 


10. Worked-out example of euclidean plane geometry 


We pose a problem in euclidean plane geometry on which to practice the theory 
developed up to now: 


Given a point A lying on an oriented line m, and a second point A' 
lying on a second oriented line m', construct the unique direct isometry 
mapping A to A' and m to m'. 


The problem is illustrated in Fig. 12 (left), including orientation on the two lines. 
We assume the points and lines are normalized, and define to begin with the 
intersection point of the lines and the joining line of the points: 


M := m A m', a := A V A 1 


The direct isometry we are seeking is either a rotation or a translation. In the 
former case, the center of rotation has to be equidistant from A and A', that is, it 
lies on the perpendicular bisector of the segment AA'. To construct this we first 



Figure 12. Left to right: the problem setting, the solution, in¬ 
terpolating the solution. 
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obtain the midpoint, and then, applying Sect. 4.4 construct the perpendicular line 
through the midpoint: 


Am :— A + A . r :— A m ■ n (— A m a) 


The condition that m maps to m' implies that the center of rotation is the same 
distance from m as from m', that is, lies on the angle bisector of the two lines. We 
choose the difference in order to respect the orientations of the lines, as the reader 
can readily confirm. The desired center is then the intersection C of r and c. 

c := m — m', C := r A c 


The final step is to construct the desired isometry. We can (for a rotation) find two 
lines through C that meet at half the desired angle of rotation: the line A V C and 
the perpendicular bisector r satisfy this condition. Then form the rotor of their 
product; the rotation is then the sandwich operator defined by this rotor. 

s := A V C, g := rs, g(X) := gXg 

One can also calculate the angle a between the two mirror lines from the equation 
cos a = r ■ s, and use this to calculate g as an exponential: g = e“ c . 

Exercise: Show that the above construction also yields valid results when C is 
ideal, and that the resulting isometry is a translation. 


11. Directions for further study 

For readers who are intrigued by the approach presented here, there are several 
natural directions for further study. If one wants to stay within plane geometry, 
there are many themes that could be cast into the PGA format. For example, one 
could explore calculus and differentiation in the plane, including point-wise and 
line-wise curves, point- and line-valued functions, etc. For a general introduction 
to differentiation in geometric algebra see |DFM07] . Ch. 8. This could lead to a 
treatment of 2D kinematics and rigid body dynamics. Or, one could use the discus¬ 
sion of three-way products in Sect.[5]as a starting point for formulating the theory 
of triangles and triangle centers in this language. One practical direction would be 
to apply the theory sketched here as a framework for 2D graphics programming. 

Another natural direction is to move from 2 to 3 dimensions and explore 
the euclidean PGA P(Rg 0 x ) for euclidean 3-space E 3 . Available resources include 
IGunlla! (Ch. 7), [G un Tlbl . and [Gun 1 Id . While many results presented here 
generalize without surprises to 3 dimensions, one conceptual challenge presented 
in moving to 3 dimensions is that the space of bivectors, crucial to kinematics 
and dynamics, is no longer exhausted by the simple bivectors (which in this case 
represent lines in 3-space); the non-simple bivectors, known classically as linear line 
complexes , exhibit much more complex - and interesting - behavior. An exhaustive 
treatment of the geometric product modeled on the one presented in the first half 
of this article would accordingly yield a richer, more complicated picture. 

Practitioners of non-euclidean geometry may be interested to know that the 
approach outlined here for the euclidean plane can be carried out analogously for 
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Figure 13. Using P(R 2;li0 ) do hyperbolic plane geometry. 


the hyperbolic and elliptic planes using the algebras P(Rg 0 0 ), resp ., ± ^ 2 i o 


Most of the features discussed above for the euclidean plane have non-euclidean 
analogies which possess a similar elegance and succinctness. An introduction to 


these metric planes is given in Ch. 6 of |Gnn 11 ai . from which Fig. 13 is taken. This 
presents a metric-neutral approach, that is, results are stated whenever possible 
without specifying the metric. 


12. Evaluation and conclusion 

We have shown that traditional euclidean plane geometry can be formulated in a 
compact and elegant form using P(M 2 o i)- We have successfully applied the algebra 
to a variety of practical problems of plane geometry and have encountered no 
obstacles to the program of extending it to all aspects of euclidean plane geometry. 

How do these results compare to existing approaches? Plane geometry is 
usually handled with a mixture of analytic geometry, linear algebra, and vector 
algebra. The foregoing has established that P(R 2 0 j) offers a variety of desirable 
“infrastructure” features which this mixed approach does not offer: 

1. It is coordinate-free (for details see Appendix |A|. 

2. Points and lines are equal citizens, rather than lines being defined in terms 
of points. 

3. Ideal elements are integrated organically, both in incidence (intersection of 
parallels) and metric relations. 

4. Join and meet operators are obtained from the Grassmann algebra. 

5. Isometries are represented by versor sandwich operators that act uniformly 
on primitives of all grades. The rotors have an exponential representation. 


12 We favor using the dual construction here also (even though it is not strictly required) since 
then reflections in lines are represented by sandwiches with 1-vectors. In the standard approach, 
where 1-vectors are points, such sandwiches represent the less familiar, less practical “point 
reflections”. 
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6. The geometric product provides a rich, interrelated family of formulas for 

distance and angle integrating seamlessly both euclidean and ideal elements. 

The last point above reflects a novel feature of P(R 2 oi) °f special note: 
euclidean and ideal elements are tightly interwoven in an organic whole. See for 
example the discussion of the 2-way products in Sect. [4] and the collection of for¬ 
mulas in Sect. [6j This tight integration is, to the best of our knowledge, available 
nowhere else. We think it deserves to be better known and understood. The dis¬ 
cussion in Sect. |A~4| below makes a modest start towards a deeper understanding. 

Implementing this algebra within modern programming languages presents 
no significant challenges. The author has implemented it in Java, JavaScript, and 
Mathematica (at different times, for different purposes) and successfully applied 
the resulting toolkit to a variety of practical geometric and graphical problems. 
The resulting infrastructure gains, in comparison to traditional approaches, have 
been gratifying. 

How does P(R* 0 -J compare to the other two geometric algebras mentioned 
at the beginning of the article? [Gal07] is a treatment of plane geometry based on 
M 2 ,o,o- While entirely appropriate as an introduction to GA at the high school level, 
it makes extensive use of non-GA techniques to overcome the limitations of R 2 ,o,o> 
which unlike the euclidean plane contains a distinguished point (the origin), and 
can by itself model neither parallelism nor translations. One of the leitmotifs of 
this article has been to show how ^ 2 , 0,0 is embedded organically within P(M 2 01 ) 
as the ideal line w, so all the features of R 2 ,o,o can be accessed easily in the model 
presented here. We are not aware of an analogous treatment of plane geometry in 
CGA to the one presented here. }Gunl6b] provides a general comparison of CGA 
and PGA for euclidean geometry and establishes that for flat geometric primitives, 
such as the domain of classical plane geometry treated in this article, PGA displays 
advantages over CGA with regard to robustness, simplicity of representation, and 
ease of learning. 

To sum up: we have demonstrated that the model of plane euclidean geom¬ 
etry provided by PGA is complete, compact, computable, and elegant. Whether 
considered pedagogically, practically, or scientifically, we believe PGA provides a 
viable alternative to traditional approaches to euclidean plane geometry. By help¬ 
ing to modernize the teaching of euclidean geometry, it could make an important 
contribution to the task mentioned at the beginning of the article, of bringing 
the dramatic advances in 19th century mathematics in geometry and algebra to a 
wider audience. 


Appendices 

Appendix A. Coordinate-free description 

We provide here a modern, coordinate-free description of the algebra instead of the 
more traditional coordinate-based approach used above in Sect. |3.2|and Sect. |3.4| 
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A.l. Foundations 

Let V be a real, 3-dimensional vector space with dual space V*. We construct 
a geometric algebra A based on V using the signature (2,0,1). We describe it 
here algebraically, and postpone until later the geometric interpretation. We begin 
by recalling some basic facts and definitions regarding the underlying Grassmann 
algebra G based on V: 

• G is a graded algebra consisting of 4 grades: 

• /\°(V) is the 1-dinrensional subspace of scalars Ml. 

■ can be identified with V. 

■ /\ 2 (V) can be identified with V*. 

■ /\ 3 (V) is a 1-dimensional vector space of pseudoscalars RI. I is defined 
more precisely below in Sect. pwPl 

• An element of /\ k (V) is called a k- vector. 

• There is an anti-symmetric bilinear product A (called the wedge or Grass¬ 
mann product) defined on G that mirrors the subspace structure of weighted 
subspaces of V. For a € f\ k (V) and b £ 

■ a A b = 0 •£=> a and b are linearly dependent. 

• Otherwise, aAb £ /\ k+m (y) represents the weighted subspace spanned 
by a and b. 

• Let a £ /\\V) and A £ /\ 2 (V). We say a and A are incident aA A = 0. 

• For a vector subspace T C f\ k (V) define the outer product null space := 
{x £ /\ 3 ~ k (V) | t A x = 0 V t £T}. 

• Notation: For a multi-vector M £ G, M = where (M)*, is the 

grade-A: part of M. 


A.2. Euclidean and ideal elements 

The inner product of the geometric algebra can be represented by a symmetric 
bilinear form B : V ® V —> R. The kernel of B is defined as: 


N := {n £ V | B( n, x) = 0 V x} 

The signature of the inner product is (2,0,1). The 1 in the third position gives the 
dimension of N. So, N is a 1-dimensional vector sub-space of V. As such, it is gen¬ 
erated by an element w, which we will specify more precisely below in Sect. A.2.2| 
Elements of N are called ideal vectors. Vectors not in N are called euclidean (or 
proper). N A consists of bivectors incident with w, and is a 2-dimensional subspace 
of A (V). An element of N A is said to be an ideal bivector; all other bivectors are 
euclidean (or proper). 


A.2.1. The square of a 1-vector; normalized euclidean vectors. In a geometric 
algebra, the geometric product is defined on 1-vectors by 

ab = ab + aAb 

where a • b = B( a, b) and a A b is the the exterior product of the underlying 
Grassmann algebra. The geometric product m 2 for a 1-vector m reduces to m • m 







26 


since the wedge product is antisymmetric. For m = w, w 2 = w- w = 0 since ui £ N. 
For any euclidean vector m, 

m 2 = m ■ m = k £ R + 

We define the norm |jmj| := y/ m 2 . Then m n := y/k m satisfies ||m„|| = 1; such 
a vector is said to be normalized. 

A.2.2. The square of a 2-vector. From the above, there are two sorts of bivectors, 
ideal and euclidean. For ideal U, U = w A m for some euclidean vector m. And, 
since w £ N, oj A m = u m. Then U 2 = —ur m 2 = 0. Using the following exercise, 
it is easy to calculate that P 2 = —1. Hence a bivector is ideal •<=>■ its square is 
zero. 

Exercise: For normalized euclidean P, one can find two orthonormal euclidean 
1-vectors m and n such that P = mn. 

A.2.3. Normalized euclidean 2-vectors. We could define a normalized euclidean 
bivector to be a bivector satisfying P 2 = —1. But we can do better, as the following 
discussion shows. Let P be any euclidean 2-vector satisfying P 2 = —1. Recall the 
definition of the scaled magnitude function S in Sect. |3.2[ We fix w to be the unique 
element of N satisfying S(u> A P) = 1, and define I := w A P. We show that these 
definitions don’t depend on P, and that the value of S(u> A P) can serve as a norm 
for bivectors. 

Lemma 1 . For euclidean bivector P and ideal bivector U, (PU)o = 0. 

Proof. Choose m £ H P^ with ||m|| = lp^Then U = Amw for A £ M*. Write 
P = nm where n is normalized and orthogonal to m. Then 

PU = (nm)(Amw) 

= An(m 2 )w 
= Ana; 

Here we have used associativity of the geometric product, and the fact that m is 
normalized. Finally, since ui £ N, (nw)o = n • w = 0. □ 

Lemma 2. Given euclidean bivectors P and Q, Q = AP + U for some A £ K, A 7 ^ 0 
and U £ . Furthermore, Q 2 = A 2 P 2 . 

Proof. The first part follows by observing that C A 2 (^0 is a subspace of co¬ 
dimension 1 in A 2 (U), and Q, P ^ N^. The second assertion follows by observing: 

Q 2 = (AP + U ) 2 

= A 2 P 2 + A(PU + UP) + U 2 
= A 2 P 2 + 2A(PU ) 0 
= A 2 P 2 
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Or define m = U V P and normalize m. 
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Here we have used the fact that the grade-0 part of the geometric product PU 
is the symmetric part of the product, that U 2 = 0 for ideal U, and the previous 
lemma. □ 

Theorem 1. Given euclidean bivectors P and Q such that P 2 = Q 2 , and lo £ N. 
Then wAP = ±w A Q. 

Proof. By Lemma 2, Q = AP + U for U € Njf. Since Q 2 = P 2 , A = ±1. Wedging 
with co yields co A Q = AP+uAU = ±co A P. □ 

The preceding theorem allows us to obtain a stronger normalization than the 
condition Q 2 = —1. Define the norm of a bivector to be ||Q|| := S(co A Q). We say 
that a euclidean bivector Q is normalized when ||Qj| = 1. In every one-dimensional 
vector subspace of f\ (V), there are two solutions {Q, Q} to Q 2 = —1. ||Q|| = 1 
picks out exactly one of these solutions. The uniqueness of this result simplifies 
many calculations. 

A.2.4. Multiplication by the pseudoscalar. Multiplication by the basis pseudoscalar 
I is an important operation, sometimes called the polarity on the metric quadric. It 
maps an element to its orthogonal complement with respect to the inner product. 
This multiplication is important enough to merit its own notation 

n(X) := X - 1 := IX 

The result is called the polar of X. By the previous section, a euclidean bivector 
Q is normalized *£=> I := wQ. Then the polar of a normalized euclidean point 
P is given by P 2 - = IP = — co since P 2 = —1. 

X 2 - is sometimes called the inner product null space of X. Note in contrast 
that Xj^ is the outer product null space. 

The situation is a little more complicated for 1-vectors. Let m be a normalized 
euclidean 1-vector. Let n be a 1-vector orthogonal to m. Then the product nm 
is a normalized euclidean 2-vector, hence I = wnm and m x := Im = wn = U, 
where U is an ideal bivector. 

Exercise: The kernel of n, restricted to 1-vectors, is N, while the kernel of n, 
restricted to 2-vectors, is N^. 

A.3. Ideal inner product on ideal bivectors 

We saw above that euclidean bivectors can be normalized, but an ideal bivector 
U satisfies U 2 = 0 hence cannot be normalized in the same way. However, there 
is a way to define an alternative norm - along with an associated inner product 
- on the ideal bivectors. We define this ideal inner product and then show how 
to derive the complete inner product structure on the euclidean elements (of all 
grades) from this ideal inner product. 

A.3.1. The quotient space V/co. Define an equivalence relation on the set of eu¬ 
clidean vectors: 

m = n <£=> 3c€K such that m — n = cco. 

Let the equivalence class of m be denoted by [m]. Define a symmetric bilinear 
form B on the resulting quotient space V/co by B([m], [n]) := B(m,n). This is 
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well-defined. For if m and n are two other representatives, then m = m + cui and 
n = n + dui. _B(m, n) = B{ m + cuj, n + du) = B( m, n) since oj £ N. 

A.3.2. An “ideal” inner product on . Furthermore, 

m = n <=$■ n(m) = n(n) 

(=>): If m e n, then m = n + clo and II(m) = II(n) + cll(a;) = II(n) since 
uj £ ker(II). {<=)'■ If n(m) = II(n), then by linearity II(m) — II(n) = I(m n) = 0. 
This means m — n £ fcer(II). Hence, by the previous paragraph, m — n = coj for 
c£l. Thus n : V/uj —» N^ defined by n([m]) := n(m) is well-defined. In fact we 
have shown that it is a bijection and hence has a well-defined inverse. 

Use this inverse to transfer the inner product B([m], [n]) onto the ideal bivec¬ 
tors via 

(U, W)oo := H(n- 1 (U),f[- 1 (W)) 

It’s not hard to show that (, ) oc is the standard positive definite inner product on 
(since we began with the signature (2,0,1)) . This induces a norm on ideal 
bivectors by ||U||oo := sj (U, U)^. It is always possible to choose a representative 
for U so that HUHoo = 1. 

A.4. Recreating the (2, 0,1) inner product from the ideal inner product 

It’s tempting to view the ideal norm (, )oo as something ad hoc added on to the 
algebra P(R5 o,i)- However, the above discussion supports the contrary interpreta¬ 
tion that the ideal inner product (,)oo on ideal bivectors is the primary structure 
from which the inner product (2, 0,1) on vectors is derived, rather than vice-versa. 
For, let a and b be two euclidean vectors, and A := a A lo and B := b A w be 
their wedge product with the ideal vector. Then define a symmetric bilinear form 
B( a,b) := (A,B)oo. From the above discussion it is clear that B is well-defined, 
and in fact, B = B. So one can begin with the ideal bivector subspace equipped 
with the signature (2,0, 0) and “push” it in this straightforward way onto the eu¬ 
clidean 1-vectors to obtain the euclidean plane. Similar constructions work for any 
dimension. 

A.5. Interpretation with respect to P(R j (| x ) 

The above treatment has been carried out for an abstract real vector space V of 
dimension 3. To arrive at the algebra P(R 2 o,i) one mu st specify V, as outlined in 
Sect. [ 2 ] above, which leads to the choice V := (M 3 )*, the dual space of K 3 . In the 
resulting vector space geometric algebra 0 n 1-vectors represents oriented planes 
through the origin and 2-vectors represent standard vectors. In the second step, 
the algebra has to be projectivized to form P (]&2 0 i)- Hence, 1-vectors transform 
to lines and 2-vectors become points. In particular, to represents a plane in (R 3 )*, 
and when projectivized represents a line, the ideal line of the euclidean plane. 
The ideal bivectors are ideal points, incident with ui. Interpreting the contents of 
Sect. |A.2.2| in this light: the difference P — Q, for normalized euclidean points P 
and Q, is an ideal point. This is reminiscent of how free vectors are defined to be 
the difference of two euclidean points. In fact, ideal points are equivalent to free 
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vectors, an insight already made by Clifford in |Cli73l . so that the vector algebra 
^ 2 , 0,0 is contained here as the ideal line with its ideal inner product. 

The equivalence classes of V/u, in the context of P(M 2 o i)> are families of 
parallel lines, which share a common ideal point. Such a set of lines is known as a 
line pencil in classical projective geometry; in this case the pencil is centered on 
(or carried by) an ideal point. To see this: m = n 4==> m — n = cut. The point 
U := m A n satisfies UAm = UAn = 0. Hence U A u> = 0, which shows that U 
is ideal, as claimed. The metric polarity n, in this context, maps an equivalence 
class [m] to an ideal point perpendicular to the ideal point U. It maps all euclidean 
points (2-vectors) to the ideal line. 

Equipped with this coordinate-free foundation of the algebra, the reader can 
now rejoin the article at Sect. [4} 
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